Switching chaos on/off in Duffing oscillator 
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Periodic forcing of nonlinear oscillators generates a rich and complex variety of behaviors, ranging 
from regular to chaotic behavior. In this work we seek to control, i.e., either suppress or generate, 
the chaotic behavior of a classical reference example in books or introductory articles, the Duffing 
oscillator. For this purpose, we propose an elegant strategy consisting of simply adjusting the 
shape of the time-dependent forcing. The efficiency of the proposed strategy is shown analytically, 
numerically and experimentally. In addition due to its simplicity and low cost such a work could 
easily be turned into an excellent teaching tool. 

PACS numbers: 05.45.Ac, 05.45.Gg, 45.80.+r, 45.40.Cc 



Almost 100 years ago, as a pioneer, G. Duffing pro- 
posed a periodic forced nonlinear oscillator [l| and with 
it began the idea that simple systems can display an ex- 
traordinary rich and complex dynamic. Today, this key 
model has a vast number of applications in science and 
engineering [2|. Forced nonlinear oscillators displaying 
chaotic behavior have been observed ubiquitously across 
a broad range ofphenomena, including fluid devices 
atomic systems [4| , cardiac tissues [5[ , chemical reactions 
For some applications, the rich and complex chaotic 
dynamics are desired (e.g., fluid mixing or data encryp- 
tion), while for others such dynamics can present major 
drawbacks and significantly decrease performance (e.g., 
permanent magnet synchronous motor Q). During the 
past half century, chaotic phenomena have been under 
intense investigation. Such booming interest has led to 
powerful tools that today enable us to elaborate control 
strategies for such rich and complex phenomena. More 
precisely, these strategies use tools such as K.A.M the- 
ory, resonance phenomena and linear stability of periodic 
orbits in order to eliminate HUTU, confine [n| (at a pre- 
scribed location and for a specific size), or spread chaos 
fl2| . In this paper, we present a facile control strategy 
that enables us to turn on/off the chaotic behavior of the 
Duffing oscillator. Our strategy consists of simply and 
only adjusting the shape of the time-dependent forcing 
to either put the system in a chaotic aperiodic or regular 
periodic oscillatory regime. 

The dynamical system considered here is the classi- 
cal Duffing's equation [1] which is one of the most com- 
mon examples in chaotic oscillation teaching texts and 
research articles. Indeed, such a dynamical system exter- 
nalizes all key characteristics of chaotic system including 
manifold intersections, strange attr actor, positive Lya- 
punov exponent, etc... (see [13| for a detailed introduc- 
tion to this system). The equation describes the motion 
of a damped oscillator with a nonlinear restoring force. 
This mechanical system models a spring pendulum with 
a nonlinear spring stiffness, the bending angle is then 
governed by: 
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FIG. 1. Two characteristic shapes of the bi-chromic forcing 
with the pairs of parameters ki — {e rn = .2,<fi = ty} and k,o — 
{e m = .9,0= .35}. 



where 5 represents the viscous damping coefficient (a, /?), 
the linear and cubic coefficient in the torsional restoring 
force. The ingenuity of this work lies in the time-periodic 
forcing 7 (t; e, cj, e m , m, <fi) which is a bi-chromic function 
of the form 



7(t;e,cj,e m ,ra, 4>) 



N(<f>) 



(cos (oji) + e m cos (moot + </>)), 



6 = a6 + p6 3 



' + 7(t;e,o;,e m ,ra,< 



(1) 



where A 7 ' is a normalizing function such that — e < 7 < e 
(see FigH]). 

The strategy consists of simply and only adjusting the 
shape of the forcing by varying a pair of parameters 
ft = {e m , <t>} i-e., the amplitude ratio and frequency while 
the frequency ratio is maintained fixed at m = 3 (see 
FigJT]) in order to turn on/off chaotic behavior. 

As a first step, it is logical to look at the unperturbed 
unforced case i.e 8 = 0, 7 = 0. Although this case is 
trivial it is an excellent groundwork to investigate the 
system studied. In this case the phase space (0, 9) is 
characterized by three, one hyperbolic unstable and two 
twin elliptic stable, equilibrium points Fig. O Almost all 
orbits are closed curves showing partial periodic oscilla- 
tions when initial energy is low (concentric orbits around 
elliptic points) or complete oscillations when initial en- 
ergy overcomes a certain threshold (orbits around the 
pair of elliptic and hyperbolic points). This threshold 
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FIG. 2. Phase space in the unperturbed case i.e, 6 — 0, 
7 = with one hyperbolic unstable (cross) and two twin 
elliptic stable (dot) equilibrium points. The characteristic 
separatrix (dotted black line), partial (solid blue line) and 
complete (solid red line) oscillations are represented. 

or separatrix corresponds to the non-closed homoclinic 
orbits that links the stable and unstable direction of the 
hyperbolic point. In our approach to see how the strategy 
works, the separatrix and more precisely the stable and 
unstable manifolds that composed it play a crucial role. 
Indeed, the ability to turn on/off chaos by adjusting the 
forcing shape (i.e., by adjusting the parameters e m and (j>) 
can be seen analytically via Melnikov approach. Indeed, 
a necessary condition for chaos is that the stable and 
the unstable manifolds of the corresponding hyperbolic 
fixed point intersect transversally. One way to verify this 
condition is to look at the simple zeros of the Melnikov 
function which corresponds to the distance between the 
stable and unstable manifolds (order of the perturbation) 
measured along a direction that is perpendicular to the 
unperturbed homoclinic orbit. In our case, the existence 
(or nonexistence) of zeros in the Melnikov function can 
adroitly be seen with the sign of 

M = max|F(£ o ;a,e,cj,e m ,ra,0)| - \R(a,f},S) \ , (2) 

where, F (to; a, e, a;, e m , m, <p) = irue/N ((f)) x 
[sech (uJ7^/2^/a) sin (uto) • • • 
• • • + me m sech (mujir /2^/a) sin (mouto + (/>)] 

is a periodic function corresponding to the forcing and, 
R(a,f3,5) = 4a 3 / 2 (5/3v // 2/3 represents the unforced 
oscillator's threshold for the possibility of chaos. In the 
case where Mo > the stable and unstable manifolds 
intersect transversally turning on the possibility of 
chaos. In the other case i.e., Mo < the manifolds do 
not intersect and the possibility of chaos is turned off. 
Figure [3] illustrates Mo vs. the forcing shape parameters 
i.e., k = {e m ,0} while all the other parameters are 
maintained fixed with values corresponding to the 
experiment Fig. [6] (i.e., a = .028, f3 = -.166, 6 = .020, 
e = .0022 and u = .20). With Fig. [3] we clearly see that 
by adjusting parameters e m and <p we are able to turn 
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FIG. 3. Mo versus the pair of parameters controlling the 
forcing shape i.e., k = {e m , </>} • 

on the possibility of chaos (from very cold to hot color) 
and vice-versa, turn off the possibility of chaos (from 
hot to very cold color). 

From these analytical results, based on Melnikov 
approach and illustrated by Mo (see Fig. [3|), two 
extreme regions in the parameter space controlling 
the forcing shape can be distinguished: one very hot 
around ^i = {e m = .2,0 = 7r}, and one very cold around 
= {^m = -9,0 = .35}. Figure [H qualitatively shows 
the ability to turn on/off chaos by numerically integrat- 
ing the phase space evolution of the system Eq. [T] for the 
pairs of parameters fti/o- 

On one hand, when chaos is turned on (i.e., parameters 
controlling the forcing shape are taken as ^i), the system 
oscillates by wandering unpredictably in a large domain 
of the phase space without any apparent pattern (see 
Fig. Ufa)). On the other hand, when chaos is turned off 
(i.e., parameters controlling the forcing shape are taken 
as fto), the system oscillates by staying trapped on a 
closed curve of the phase space with a periodic pattern 
(see Fig.EIJb)). 

In order to determine if a system is chaotic and quan- 
tify its chaoticity it is important to study the evolution of 
a cloud of neighboring initial conditions. More precisely, 
we want to quantify the sensitivity to initial conditions 
which is a key ingredient in chaotic systems. This sensi- 
tivity to initial conditions can be quantified as: 

\6X(t)\*e M \5X(0)\, 

where (t) is the evolution of a close cloud of neighbor- 
ing initial conditions and A is the mean rate of separation 
of the system called Lyapunov exponent. When A > 0, 
any two orbits that start out very close to each other 
diverge quickly, otherwise when A < all points in the 
neighborhood converge toward the same periodic orbit. 
In order to compute the Lyapunov exponent we first con- 
sider the time evolution of the Jacobian J 1 ( 9 ) given 
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FIG. 4. Phase space from numerical simulation when chaos is 
turned on (i.e., e m = .2, (f> = tt) and off (i.e., e m = .9, <j) — .35) 
(a)-(b), respectively. 



by the tangent flow and the matrix of variations as 



~dt 



1 

dO d Q 



where J° = / is the two-dimensional identity matrix. 
Then, the leading Lyapunov exponent for a finite-time 
t = r is defined as A (r) = ^ In |A ma:E | where A max is the 
largest eigenvalue (in norm) of the Jacobian J T . 
Figure [5] shows the leading Lyapunov exponent A vs. time 
t when chaos is turned on/off (i.e., with forcing shape 
parameters fti/o) • On one hand, when chaos is turned 
on Fig. [5] (a) displays a Lyapunov exponent converging 
toward a positive value of A = .015. Such a value indi- 
cates high sensitivity to initial conditions in the system 
and consequently high chaotic behavior. In the other 
hand, when chaos is turned off, Fig. [5] (b) shows a nega- 
tive Lyapunov exponent value indicating no sensitivity to 
initial conditions and consequently no chaotic behavior. 
All orbits converge toward the periodic orbits shown in 
Fig-IUJa). After having been applied analytically and nu- 
merically, this strategy for turning on/off chaos has also 
been implemented experimentally. The design and con- 



struction of our experimental system based on [14| and 
is illustrated in Fig. [6l Such a setup presents the advan- 
tage of being economical and simple to realize while still 
being able to be qualitatively described by Eq. (pQ). The 
oscillator consists of a vertical steel beam with a brass 
weight clamped at its top. This beam is subjected to a 
forcing produced by a pair of electromagnets connected 
to an arbitrary wave generator. This generator fed the 
electromagnets with a bi-chromic voltage wave identical 
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FIG. 5. Leading Lyapunov exponent A vs time t from numer- 
ical simulation when is turned on (i.e., e m = .2, <j> = tt) and 
off (i.e., Cm = -9, <j> = .35) (a)-(b), respectively. 
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FIG. 6. Sketch of experimental setup. 

(a) Mechanical oscillator with weight, steel beam [l5j], strain 
gauges [16[] and electromagnet s |17ll from top to bottom. 

(b) Arbitrary waive generator [18| . 

(c) Sketch of the differential amplifier electrical [191 ]. 

(d) Data acquisition [20[| and signal treatment system. 



to 7 (see Fig.[Tj). In addition, four strain gauges glued on 
the beam send an electrical signal to an differential am- 
plifier electrical circuit providing voltages proportional to 
the bending angle from the vertical and the angular ve- 
locity 0. The experimental procedure consists of starting 
the oscillation around (0q = 0, #o = 0), then waiting 20 
min (i.e., more than 200 times the forcing period) in order 
to pass the transient mode to make sure the system has 
well converged toward its attr actor, and finally recording 
the data at 60 Sa.s -1 for 20 min. Figure shows exper- 
imental results for the phase space evolution of our ex- 
perimental Duffing oscillator Fig. [6j These results verify 
the strategy proposed i.e., by simply and only adjusting 



4 



(a) 
.03 



.03 




Due to the universality of the Duffing equation in sci- 
ence, we expect this strategy to be applicable to other 
oscillatory experimental systems as well. 

The authors are grateful to C. Chandre for his help 
and suggestions, and thanks M. Labour, N. Zaharia, B. 
Belda and H. Wang for their assistance in performing the 
experiments. 
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FIG. 7. Nondimensionalized phase space from experiment 
when chaos is turned on (i.e., e m = .2, <fi = tt) and off (i.e., 
e m = .9, <j> = .35) (a)-(b), respectively. The nondimension- 
alized experimental parameter values are the same as in Fig. 
31 These nondimensionalized values correspond to: a forcing 
frequency of 18 mHz and an amplitude of 1.25 Vpp on the 
magnets; a torsional restoring force cubic and linear coeffi- 
cient of 2.5 and 15.1 mN.m, and a torsional viscous damping 
coefficient of .31 mN.m.s. 



the forcing shape we are able to turn on/off chaos. As 
predicted by the analytical and numerical results, chaos 
is turned on/off when the forcing shape parameters are 
adjusted as ^!/ (see Fig. [T] for illustration). On one 
hand, when chaos is turned on, the system oscillates in 
unpredictable fashion in a large part of the phase space 
domain. On the other hand, when chaos is turned off, 
the system oscillates in a completely predictable peri- 
odic pattern. We have shown analytically, numerically 
and experimentally that chaotic behavior of the famous 
Duffing oscillator can be controlled successfully by the 
bi-chromic forcing 7. More precisely, by adjusting the 
forcing shape (i.e., the amplitude ratio e m and the phase 
difference qb between the two harmonics composing the 
forcing) chaotic behavior can be turned on/off. This con- 
trol strategy has the great advantage of being extremely 
easy to set up since neither feed back loop nor addition 
of time variation input are necessary. Furthermore, this 
strategy does not require important investments such as 
increase in force amplitude or frequency. 
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